Strong Color Field Baryonic Remnants in Nucleus-Nucleus Collisions at 

200A GeV. 
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The effects of strong color electric fields (SCF) on the baryon production at RH1C are studied 
in the framework of HIJING/BB (v2.0) model. The particle species dependence of nuclear modi- 
fication factors (NMF) are analyzed for Au+Au collisions at 200A GeV. A doubling of the string 
tension leading to a modification of the strangeness suppression according to Schwinger mechanism 
is shown to provide an alternate explanation to coalescence models for the interpretation of the 
observed baryon and meson production at moderate pr and results in a predicted enhancement in 
the (multi)strange (anti)hyperon production. 



PACS numbers: 25.75.Dw; 24.10.Lx. 
I. INTRODUCTION 

While the phase transition from hadronic degree of 
freedom to partonic degrees of freedom (quarks and glu- 
ons) in ultra-relativistic nuclear collisions is a central fo- 
cus of recent experiments at the Relativistic Heavy Ion 
Collider (RHIC), data on baryon and hyperon produc- 
tion has revealed interesting and unexpected features at 
RHIC that may be of novel dynamical origin instead. 
The so called baryon/meson anomaly [1-4] is observed as 
a large enhancement of the baryon to meson ratio and a 
large difference of the nuclear modification factor (NMF) 
between total charged [5] and neutral pions (ir°) [6] at 
moderate transverse momenta (2 < pr < 5 GeV/c). 

In a previous paper [7] we studied the possible role of 
topological baryon junctions [8] in nucleus- nucleus colli- 
sions. We have shown in the framework of HIJING/BB 
v2.0 model, that junction-antijunction (J J) loops with an 
enhanced intrinsic transverse momentum kx ~ 1 GeV/c, 
a default string tension kq = 1 GeV/fm, and a di-quark 
suppression factor (PARJ(l)=7 gg =0.07) provide a par- 
tial explanation of the baryon/meson anomaly [7]. That 
model therefore provides an alternative dynamical expla- 
nation of the data to recombination models [9]. Within 
HIJING/BB v2.0 [7] one of the main assumptions is: the 
strings could survive and fragment [12,13], and in partic- 
ular populate the mid to low pr range. In contrast, in 
the recombination picture [9] or in hydrodynamical ap- 
proach [14] all coherent strings are assumed to become 
rapidly incoherent resulting in rapid thermalization. 

In this paper we explore further dynamical effects as- 
sociated with long range coherent fields (i.e strong color 
fields, SCF) including baryon junctions [8] and loops 
[10] that may arise in nuclear reactions. Our empha- 
sis here will be on the novel baryon observables mea- 
sured at RHIC. In nucleus-nucleus collisions the color 
charge excitations may be considerably greater than in 
nucleon-nucleon collisions due to the almost simultane- 
ous interaction of several participating nucleons in a row 



[11,12]. Molecular dynamics models [15-17] have been 
used to study the effects of color ropes as an effective 
description of the non-perturbative, soft gluonic part of 
QCD [18-20]. Strangeness enhancement [21-29], strong 
baryon transport [30], and increase of intrinsic kx [19] 
are all expected consequences of SCF. This can be mod- 
eled in microscopic models as an increase of the effective 
string tension that controls the qq and qqqq pair creation 
rates and strangeness suppression factors [11]. 

For a uniform chromoelectric flux tube with field (E) 
the probability to create a pair of quarks with mass (m), 
effective charge (e), and transverse momentum (pt) per 
unit time per unit volume is given by [31] : 
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The integrated probability {P m ) reproduces the classical 
Schwinger results [32] , derived in spinor quantum elec- 
trodynamics (QED) for e + e~ production rate, when the 
leading term in Eq. 2 is taken into account, i. e.: 
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In general in microscopic string models the heavier flavors 
(and di-quark) are suppressed according to Schwinger 
formula [32]: 
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where K — \eE\ is the string tension; ttlq is the effec- 
tive quark mass; (Q=s for strange quark; Q=qq for a 
di-quark), and q=u,d are the light nonstrange quarks. 

In the case of quark-gluon plasma (QGP) creation it is 
necessary to modify the dynamics of particle vacuum pro- 
duction at short time scales, and the abundance of newly 
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produced particle may deviate considerably from the val- 
ues obtained for the constant field [33] . Two possible pro- 
cesses may lead to an increase of strangeness production 
within the framework of Schwingcr mechanism arc: i) in- 
creasing the field strength by a modified string tension 
(k) [11,19,25,26,34], or ii) dropping the quark masses due 
to chiral symmetry restoration [35,34,36,37]. A specific 
chiral symmetry restoration could be induced by a rapid 
deceleration of the colliding nuclei [38]. 

Present estimates [39] of the current quark masses 
range from: m u =1.5-5 MeV; m d =3-9 MeV, m s =80- 
190 MeV. For di-quark we consider m gi j=450 MeV [40]. 
Taking for constituent quark masses of light non-strange 
quark M u4 = 230 MeV, strange quark M s =350 MeV [41], 
and di-quark mass M qq =550 ± 50 MeV as in Ref. [40], 
it is obvious that the masses of (di) quark and strange 
quark will be substantially reduced at the chiral phase 
transition. If the QGP is a chirally restored phase of 
strongly interacting matter, in this picture the produc- 
tion of strange hadrons will be enhanced [37]. In this 
case a possible decrease of the strange quark mass would 
lead to a similar enhancement of the suppression factors, 
obtained (in microscopic models) by an increase of string 
tension [18,19,34,33,36]. Moreover, if we consider that 
Schwinger tunneling could explain the thermal character 
of hadron spectra and that, due to SCF effects the string 
tension value k fluctuates w e can define an apparent tem- 
perature T = ^< k > /2ir [42] . 

II. OUTLINE OF THE HIJING/BB V2.0 MODEL 

Our analysis are performed in the framework of HI- 
JING/BB v2.0 model [7] that is based on HIJING/BB 
vl.10 [43]. Multiple hard and soft interactions proceed 
as in HIJING vl.37 [44]. In HIJING/BB v2.0 we intro- 
duced [7] the possible topology with two junctions [47], 
and a new algorithm where J J loops arc modeled by an 
enhancing di-quark p T kick characterized by a gaussian 
width of a qq ' = f ■ a qq = 1.08 GeV/c, with a qq — 0.360 
GeV/c (consistent with PYTHIA [45] default value), i.e. 
/=3, which we fit to best reproduce the observed pr 
spectrum of the baryons. 

Following the equations above, we take into account 
SCF in our model by an in medium effective string ten- 
sion k > k , which lead to new values for the suppres- 
sion factors, as well as the new effective intrinsic trans- 
verse momentum Ut [18-20]. This includes: i) the ratio 
of production rates of di-quark to quark pairs (di-quark 
suppression factor), ^ qq = P(qqqq)/ P(qq), n ) the ratio 
of production rates of strange to nonstrange quark pairs 
(strangeness suppression factor), j s = P(ss)/ P(qq), iii) 
the extra suppression associated with a diquark contain- 
ing a strange quark compared to the normal suppression 
of strange quark (7^), j us = (P(uuss)/P(uudd))/(j s ), 
iv) the suppression of spin 1 diquarks relative to spin 
ones (appart from the factor of 3 enhancement of the for- 



mer based on counting the number of spin states), 710, 
and v) the (anti)quark (a q = ^Jk/kq ■ a q ) and (anti)di- 

quark (a qq = a/ k/kq ■ f ■ <7 qq ) gaussian width. These 
parameters correspond to 7 g9 =PARJ(l), 7 S =PARJ(2), 
7 US =PARJ(3), 7i =PARJ(4), and <r qq = PARJ(21) of the 
JETSET7.3 subroutines [45]. Our calculations are based 
on the assumption that the effective enhanced string ten- 
sion (k), in both basic ropes {q n — q n ) and junction ropes 
(q n — q n — q n ) are the same. For elementary n strings 
and junctions this ansatz is supported by baryon studies 
[46] . A different approach to baryon production without 
baryon junctions has been proposed in [25] where SCF 
from string fusion process can lead to (qq)e — (qq)e with 
about double the string tension. Both types of SCF con- 
figurations may arise but predict different rapidity depen- 
dence of the valence baryons. We consider in this version 
2.0 of HIJING BB only baryon junction rope loops. 

There is a debate in the study of qqq system on the 
shape of A-like geometry and V-like geometry [47], [48], 
and on the stability of these configurations for the color 
electric fields [49] . In both topologies we expect a higher 
string tension than in an ordinary qq string (ky = v3 Ko 
and ka = (3/2) Ko). It was shown [49] that the total 
string tension has neither the Y nor the A-likc value, 
but lies rather in-between the two pictures. However, the 
Y configuration appears to be a better representation of 
the baryons . If two of these quarks stay close together, 
they behave as a di-quark [48] . In a dual superconductor 
models of color confinement for the Y-geometry the flux 
tubes converge first toward the centre of the triangle and 
there is also another component which run in opposite 
direction. They attract each other and this lower the 
energy of Y-configuration [47] . 

Phcnomcnological applications are currently based on 
Regge trajectory which gives the appropiate relationship 
between the mass M of the hadrons and its spin J s : J s — 
a + a s M 2 , where a ~ 0.5 is the Regge intercept; a s 
is the Regge slope. The value of the Regge slope for 
baryons is a s ~ 1 GeV~ 2 [43] that yields a string tension 
( related to the Regge slope, kq — l/2na s [50]) kq «1 
GeV/fm. This value is taken in our calculations within 
HIJING/BB v2.0 toghether with the following values for 
the suppression factors correponding to the constituent 
quark masses given above: j qq =0.02; 7 S =0.30; 7„ s =0.40; 
7io=0.05. A broadening in kx is chosen as: ct 9 =0.360 
GeV/c and o' qq = 1.08 GeV/c (i.e f=3). HIJING/BB 
v2.0 predictions using this set of parameters are labeled 
here by w/o SCF, i.e. without strong color field, or by 
k =1 GeV/fm. 

The multi-gluon exchange processes dominated by 
Pomeron exchange in high energy nucleus-nucleus col- 
lisions could be described by a Regge trajectory with 
a smaller slope a s 1 ~ 0.45 GeV -2 [51], leading to an 
increase of string tension to k k 2 k$ [19], correspond- 
ing to an increasing values for the suppression factors: 
7^=0.12; 7 S '=0.55; 7„ s '=0.63, 710 '=0.12 as well as of 
a broadening in kr , (7^=0.500 GeV/c and a qq = 1.5 
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GeV/c. The results obtained with this set of parameters 
are labeled here w/SCF, i.e. with strong color field, or 
by k= 2 GeV/fm. We note, that the increase of "intrin- 
sic fcy" (the gaussian width a'^ q ), is strongly supported 
by recent experimental values reported by PHENIX [52], 
which show an increase from p+p (fc^ J/ =1.08 ± 0.05 
GeV/c) to d+Au collisions 0^=1.36 ± 0.07 ± 0.12 
GeV/c) at ^5^=200 GeV. 



III. NUCLEAR MODIFICATION FACTORS 

Here we will concentrate our analysis on species de- 
pendence of the nuclear modification factors (NMF) Raa 
and R cp in Au+Au collisions at a/s_/vat=200 GeV. Raa 
is the ratio of the heavy-ion yield to the pp cross sec- 
tion normalized by the number of binary collisions, while 
R cp is the ratio of scaled central to peripheral particle 
yield and are defined as in Ref. [7] . For p+p interactions 
we have used in our calculations, the set of optimized 
parameters from PYTHIA [45] or HIJING [44]. 

Figure 1 shows the predicted NMF Raa for the sum 
of hadrons for two centralities. The results are compared 
to the data obtained by the PHENIX collaboration [53] . 
The data at both centralities could not be described by 
assuming only a broadening of the intrinsic kx from its 
standard value (dotted histograms) to (7^=1.08 GeV/c 
(i.e 7=3, dashed histograms). The introduction of SCF 
has an effect on the predicted nuclear modification fac- 
tors of the total inclusive hadrons and results in a better 
agreement with data (solid histograms). The data in- 
dicate at most a small variation with centrality of the 
factor "/" consistent with the broadening originating at 
the parton level. 
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FIG. 1. (Color online) Comparison of HIJING/BB v2.0 
predictions for Raa of total inclusive charged hadrons, in 
central (0-10% -left) and peripheral (60-90%-right) Au+Au 
collisions at y / sjviv=200 GeV. The results are with (solid his- 
tograms) and without SCF (dashed histograms). The label 
f3 stands for model calculations assuming /=3. Dashdotted 
histograms are the results without SCF and f=l (label fl). 
The data are from PHENIX [53]. Only statistical error bars 
are shown. 

In order to better quantify possible effects of strong 



color field on particle production we investigate species 
dependence of NMFs Raa(pt) f° r central collisons 
(Fig. 2a-d.), where higher sensitivity to SCF is expected. 
Because of their dominance, the production of pions is 
only moderately modified when we consider an increase 
of the string tension value, since the total energy is 
conserved. Taking into account SCF effects (solid his- 
tograms) results in changes at moderate px of less than 
~ 20% for the pion yield (Fig. 2a). The scaling behaviour 
in Raa(pt) of the pions is different from those of the sum 
of protons and anti- protons (Fig. 2c). The pions yield in 
central events is strongly suppressed compared to binary 
collisions scaling (Raa{pt)=^)- 
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FIG. 2. (Color online) HIJING/BB v2.0 predictions for 
species dependence of NMF (-Raa) in central (0-10%) Au+Au 
collisions at y / ijvjv=200 GeV: (a) for charged pions, (b) kaons, 
(c) inclusive p+p, (d) inclusive A + A. The results are with 
(solid histograms) and without SCF (dashed histograms). 
The dotted histograms are the predictions assuming «j=1.5 
GeV/fm. The data are from PHENIX [54] and STAR [55,56] 
collaborations. Only statistical error bars are shown. 

The hadron production in HIJING/BB v2.0 is mainly 
from the fragmentation of energetic partons. Thus, the 
observed suppression of pions in central collisions may 
be a signature of the energy loss of partons during their 
propagation through the hot and dense matter (possi- 
bly QGP) created in the collisions, i.e jet quenching. On 
the contrary, strange particles are highly sensitive to the 
presence of SCF. Kaons (Fig. 2b) and lambda (Fig. 2d), 
show an increase by a factor of two and ten, respectively. 
Such an increase results in a predicted enhancement of 
the lambda yield relative to scaled binary collisions as op- 
posed to the strong suppression predicted (and observed) 
for pions. These results are consistent with the PHENIX 
data [54] and preliminary STAR data [55,56]. In or- 
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der to study the sensitivity to string tension values we 
also present the results corresponding to an intermediate 
value of the string tension i.e. «;= 1.5 GeV/fm (dot- 
ted histograms). The discrepancy seen at low pt < 1.2 
GeV/c (solid histograms) comes from a sizeable contribu- 
tion from radial flow, not included in our model. The low 
pr region at RHIC seems to be better described within 
hydodynamical approach [14]. 

The predicted particle dependence is due to the inter- 
play between nuclear effects such as jet quenching and 
shadowing and fluctuations of the chromoelectric field in 
the early phase of the reaction. Insight into the source 
of the particle dependence is obtained from Fig. 3 and 
Fig. 4. The comparison of the NMF R cp (pt) for protons 
and pions (Fig. 3a) as well as for kaons (Fig. 3c) presents 
a behavior similar to Raa{pt), i.e shows a meson/baryon 
anomaly that is well described and was interpreted in [7] 
as due to a possible exotic gluonic mechanism (JJ loops). 
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FIG. 3. (Color online) HIJING/BB v2.0 w/SCF (k = 2 
GeV/fm) predictions for species dependence of NMFs R cp (pt) 
in Au+Au collisions at v / sjvjv=200 GeV. The results are for 
scaled (0-10%)/(60-90%). The data are from PHENIX (filled 
symbols) and from STAR (open symbols). PHENIX data 
are scaled (0-10%)/ (60-90%) [2]. STAR data are at slightly 
different centralities: scaled (0-5%)/(60-80%) for kaons and 
A's [57], and scaled (0-5%)/(40-60%) for 5 (preliminary data) 
[55]. Only statistical error bars are shown. 

In Ref. [9] it is suggested that the behavior of R cp (pt) 
may be interpreted as due to the competition between 
recombination and parton fragmentation. The results 
obtained for R cp {pt) for the strange (Fig. 3b) and multi- 
strange particles (Fig. 3d) show a small suppression rel- 
ative to binary scaling, consistent with the experimental 
results for A [2,57], and preliminary data for 5 [55]. 

In contrast, the equivalent predictions of Raa(pt) for 



protons and (multi)strange particles (Fig. 4) show a pre- 
dicted strong enhancement of NMF due to SCF effects 
dependent on the mass and strangeness content of the 
produced particles. 
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FIG. 4. (Color online) HIJING/BB v2.0 predictions with 
SCF (right) and without SCF (left) for species dependence of 
Raa(pt) in central (0-5 %) Au+Au collisions at V / Sjvjv=200 
GeV. The preliminary STAR data are from [55,56]. Only 
statistical error bars are shown. 

In particular, the model predicts a dramatic enhance- 
ment in the multi-strange (anti-)hyperon production at 
moderate pt, up to a factor of roughly 2 (relative to bi- 
nary scaling) for Ss and up to a factor of more than 10 for 
fis (Fig. 4b). The striking difference between Raa and 
R cp could be explained in our model as a consequence of 
SCF which manifest in both central and peripheral colli- 
sions. This shows that there is a clear difference between 
using peripheral Au+Au yields (as in R cp ), or p+p yields 
(as in Raa) as base-line for comparison with binary col- 
lisions scaling (Raa=Rcp=^)- 



IV. SUMMARY AND CONCLUSIONS 

In summary, we studied the influence of possible strong 
longitudinal color fields in particle production in heavy- 
ion collisions. We modeled SCF effects within HI- 
JING/BB v2.0, by varying the effective string tension 
that controls the qq and qqqq pair creation rates and 
strangeness suppression factors. We show that junction- 
(anti)junction loops and a higher string tension k = 2 n 
(k sa 1 GeV/fm), could be important dynamical mech- 
anisms in the solution of the observed baryon/meson 
anomaly. Our approach has the advantage of correlating 
many observables in the same dynamical model includ- 
ing pions, kaon, and baryons productions and of predict- 
ing species dependence of the nuclear modification fac- 
tors Raa(j?t) an d R cp {pt)i and the transverse momenta 
over a large px range. A greater sensitivity to SCF ef- 
fects is predicted for the nuclear modification factors of 
(multi)strange hyperons. In particular, the measurement 
of fl and yields would provide an important test of the 
consistency of SCF and baryon junction mechanisms at 
RHIC. 
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